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INTRODUCTION 

Many factors influence the amount of gas required to pressurize a cryo- 
genic propellant tank during the period of outflow. Besides the tank volume 
and the temperature and pressure of the incoming gas, other factors such as out- 
let flow rate, gas-to-wall heat-transfer coefficient, mass and specific heat of 
the tank wall, and the gas specific heat must be considered. A systematic ex- 
perimental investigation of these individual factors is very difficult for liquid 
hydrogen. It is desirable, therefore, to attempt analytically to determine the 
relative significance of the various parameters. 

An analysis of the tank pressurization problem for a cylindrical tank was 
made at Lewis Research Center (ref. l). A simple one -dimensional model was used, 
based on a rather restrictive set of physical assumptions. Even for the simple 
model the resulting differential equations were quite complex and a numerical 
solution was clearly indicated. The details of the numerical solution were worked 
out and a computer program was developed. Results of the analysis were compared 
with experimental results for a number of cases and the agreement was shown to 
be surprisingly good in view of the ' restrictive assumptions. 

The good agreement appeared to justify the use of the computer program for 
investigating systematically the various parameters affecting the pressurization 
problem. This investigation was carried out and the results are presented in de- 
tail in a forthcoming report (ref. 2). A brief discussion of these results and 
of the assumptions involved is given in the present paper. 


ANALYSIS 

The analysis is restricted to the cylindrical portion of the tank (fig. l) 
and only the period of time during which outflow occurs is considered. Cer- 
tain assumptions are made in an attempt to simplify the analysis and shorten the 
subsequent numerical solution while still retaining the most important features of 
the problem. A list of the assumptions and a discussion of their validity follows: 

(1) The ullage gas is nonviscous. 

(2) The velocity of the ullage gas is parallel to the tank axis and varies 
only in the axial direction. 
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(3) The tank pressure varies only with time. 

(4) The ullage gas temperature varies only in the axial direction. 

(5) The tank wall temperature varies only in the axial direction. 

(6) No heat is transferred axially in either the gas or the wall. 

(7) No condensation or evaporation occurs. 

(8) No heat is transferred at the liquid interface or at the top of 
the tank. 

With these assumptions the problem is reduced to a one -dimensional, nonsteady, 
nonviscous flow of the ullage gas with heat transfer to the tank wall. 

Although the problem is clearly not one -dimensional (radial flow must 
take place as the gas enters the tank), it is necessary to simplify the 
equations. Therefore, assumptions (l) and (2) stipulate that the pressurizing 
gas enters the tank uniformly at x = 0 (fig. l) and proceeds downward with a 
velocity that varies with time and axial location only; that is, no mixing of 
the ullage gas occurs. 

Assumption (3) is likely to be satisfied closely because of the low gas 
density and small change in gas momentum from top to bottom of the tank. 

Assumption (4) arises from experimental results obtained at Lewis with a 
cylindrical tank having a low heat leak. The assumption may not be valid for 
other circumstances. 

Assumption (5) is adequate for thin metal tank walls. 

Assumption (6) arises from the low conductivity of the ullage gas and 
the small thickness of the tank wall. 

Assumption (7) appeared to be justified by early data taken at Lewis. 
Recently taken data, however, put the assumption in doubt.. More experimental 
results, especially on larger tanks, are needed to evaluate this assumption 
properly. 

Assumption (8) has not been verified. There are likely to be some cases in 
which the heat transfer to the top of the tank, at least, cannot be ignored. 

With these assumptions, the differential equations that govern the pressuri- 
zation problem can be written (see ref. 1 for details) 
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(All symbols are defined in the appendix.) In addition to these three differen- 
tial equations in the three unknowns T, T w , and u, the following initial and 
boundary conditions are also required to determine a solution: 


(1) At the start of outflow, the gas and wall temperature distributions 
must be given. 

(2) The variation during outflow of the incoming gas temperature, the tank; 
pressure, the outlet flow rate, and the gas and wall temperatures at the inter- 
face must be prescribed. 


Furthermore, the heat-transfer coefficient must be supplied, either by an 
equation relating it to fluid properties or by using appropriate experimental 
values . 


NUMERICAL SOLUTION 

A finite difference solution of equations (l) to (3) was programmed in 
Fortran IV for use on an IBM 7094-11 computer. Backward difference equations 
were used resulting in a nonlinear set of algabraic equations that were explicit 
in the unknown variables. 

The time step At is related to the space step Ax by the requirement 

that 

Ax 

At = — rry 

where uq/t) is the velocity of the liquid surface. This restriction on At is 
used to keep the net spacing Ax constant as the solution progresses. (it is 
not a condition for stability of the numerical solution and it does not result in 
unusually small values of At). The program has been run over a very wide range 
of problems and no numerical instability has been encountered. 

The output of the computer program is the distribution of gas and wall tem- 
peratures at any desired time during outflow. The pressurant mass required at 
each instant is also determined. A typical solution uses about 200 net points in 
the x-direction for covering the entire length of the tank. The 19 solutions pre- 
sented in reference 1 averaged 24 seconds of computer time per solution. 
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EXPERIMENTAL AND CALCULATED RESULTS 

In reference 3 the authors report some of the results of a systematic 
series of liquid hydrogen expulsion experiments. The tank used was 27 inches 
in diameter and 89 inches in overall length with dished head ends. A gas dif- 
fuser was used at the inlet. The tank was constructed of 5/l6-inch 304 stainless- 
steel plate and was vacuum jacketed. The instrumentation, described in detail 
in reference 3, provided a relatively significant heat sink in some of the experi- 
ment s . 

Ten experiments (some of which were not discussed in ref. 3) were selected 
to check the analysis. These covered a wide range of outlet flow rates, tank 
pressures, and inlet gas temperature variations. Helium was used to pressurize 
in four of the cases. The detailed input data necessary to carry out the cal- 
culation is given in reference 1 for each of the experiments. Some of the prin- 
cipal data are given in table I. 

TABLE I; - LEWIS EXPERIMENTAL DATA 

^ Example, Pressure, (Outflow time, j Outflow rate,', Experimental average' Pres sur- 
/ . ! ! / 



lb/sq in. ' 

i 

sec 

■ cu ft/sec 

heat -transfer 
coefficient, 
Btu/(sq ft)(hr)(°R) 

izing 

gas 

1 

160 

350 

0.0669 

13.75 

h 2 

2 

161 

93 

.2375 

12.25 

h 2 

3 

57 

284 

.0780 

7.09 

h 2 

4 

58 

101 

.2238 

6.67 

h 2 

5 

164 

95 

.2340 

11.34 

h 2 

6 

40 

88 

0.2550 

5.13 

h 2 

7 

159 

355 

.0634 

12.31 

He 

8 

159 

90 

.2598 

11.15 

He 

9 

159 

100 

.2365 

10.45 

He 

10 

40 

309 

.0703 

5.25 

He 


Figure 2 shows the gas and wall-temperature distributions calculated at the 
end of outflow and the corresponding experimental values for each example. The 
agreement generally is good. 

Reference 4 reports the results of hydrogen experiments carried out at 
Lockheed-Georgia Company using a 40 -inch-diameter test tank 100 inches in overall 
length. The test tank was 0.090 -inch-thick stainless steel and was enclosed in a 
60-inch-diameter vacuum-tight carbon steel tank. A gas diffuser was in the top 
and an antivortex baffle was in the bottom. Perforated conical slosh baffles were 
located at various axial positions. The heat sink effect of the internal hard- 
ware could not be well estimated from the information reported. 

Nine tests are reported in reference 4 for which the system vacuum was main- 
tained. These cover two values of inlet gas temperature and a range of values of 
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initial ullage . The outflow time and tank pressure varied only slightly from 
test to test. Helium was used to pressurize in one case. Sloshing of the liquid 
was induced in all but one case. The detailed input data for the calculations is 
given in reference 1. Some of the principal data are shown in table II. 


TABLE II. - LOCKHEED -GEORGIA EXPERIMENTAL DATA 


Example 


Pressure, 


t n " 

Outflow time, | Outflow rate, ( Experimental average. Pressur- 
izing 
gas 


! lb/ sq in. 

1 sec 

cu ft/sec 

heat-transfer 
coefficient, 
Btu/(sq ft)(hr)(°R) 

1 

45.5 

89 

a 0 .672 

11.5 

2 

47.6 

103 

.560 

b 12.0 

3 

46.5 

120 

.511 

11.3 

4 

46.5 

87 

.607 

, 12.0 

5 

45.5 

99 

.609 

12.1 

6 

1 47.0 

95 

.644 

12.3 

7 

45.0 

111 

.530 

11.8 

8 

46.2 

97 

.632 

11.7 

9 

45.5 

105 

.565 

l ■ - 

13.9 

i .. 

a Flow rates are computed from reported outflow time, tank volume, 


H 2 

h 2 

% 

h 2 

He 

h 2 

h 2 

H 2 

Ho 


initial ullage. 

■'Estimated value; not given in reference 4. 


For the Lewis and the Lockheed-Georgia experiments pressurant mass require- 
ments were obtained from the analysis. Table III shows these calculated values 
along with the experimental value in each case. The percent difference is 
also shown. The average difference for the Lewis experiments is about 5 percent. 
The average difference for the LockheedKJeorgia experiments is about 4— percent. 
This agreement is better than might be expected from the simple description of 
the problem used for the analysis. 


TABLE III. - PRESSURANT MASS REQUIREMENTS 


Example (Experimental [Calculated | Percent* Example [Experimental • Calculated I Percent 



mas s , 

mass, 

j differ- 


: mass. 

mass, 

differ- 


lb 

lb 

ence 


lb 

lb 

ence 


Lewis 

data 



Lockheed-Georgia data 


1 

3.98 

i 3.95 

I -0.75 

1 

2.61 

2.81 

7.67 

2 

2.72 

2.60 

-4.41 

2 

2.13 

2.24 

5.17 

3 

1.76 

; 1.68 

: -4.54 

3 

2.86 

3.05 

6.64 

4 

1.24 

1.27 

! 2.42 ' 

4 - 

| 2.57 

2.65 

3.11 

5 

3.76 

3.51 

-6.65 

5 

5.79 

5.89 

1.73 

6 

.83 

.93 

' 12.04 , 

6 

2.47 

2.58 

4.45 

7 

8.14 

7.61 

-6.51 : 

7 

2.81 

2.86 

1.78 

8 

5.59 

5.57 

-.36 : 

8 

2.81 

2.95 

4.98 

9 

9.24 

8.48 

-8.23 

9 

2.88 

3.00 

4.17 

10 

2.70 

2.56 

-5.18 ,j 


. ] . _ ... .. 



1 

1 
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It should he noted that experimental average values of heat-transfer coeffi- 
cient were used and that the gas and wall temperature distributions at the start 
of outflow were obtained from the data. The variation of inlet gas temperature 
with time at the position x = 0 is also from the experiments. 


PARAMETRIC ANALYSIS 

The agreement shown between calculated and experimental values of pressurant 
mass requirement in the preceeding section encourages the use of the analytical 
method for examining the effect of the various parameters entering the pressuri- 
zation problem. A method for doing this and the results obtained are described 
briefly in this section. 


Dimensionless Parameters 

The following additional assumptions are made to simplify the differential 
equations (l) to (3) and the initial and boundary conditions: 

(9) The ullage gas is a perfect gas with constant specific heat. 

(10) The gas -to -wall heat-transfer coefficient is constant in space and time 
for a given example . 

(11) The inlet gas temperature, the tank pressure, and the outflow rate are 
constant . 

(12 ) The gas and wall temperatures at the liquid interface are constant and 
equal throughout the outflow period. 

(13) The gas and wall temperatures at the start of outflow are equal and 
vary linearly in the direction of the tank axis from the temperature at the 
liquid interface to a temperature at the top of the tank equal to the average 
of the inlet gas temperature and the liquid surface temperature. 

The last assumption is obviously an arbitrary choice for the initial gas and 
wall temperatures. The effect of this assumption and the others will be consid- 
ered later. 

Using these assumptions and introducing the dimensionless variables 



/s u 
U “ U L 
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a m 

T = - 7 - 
T g 


Tv = 




into equations (l) to (3) gives 


DT , A /N 

^ = 2 St_(T - T)T 
Dt g w ' 


(4) 


(5) 


where 


WX W ,/s /s 

= S W “ T w ) 
du 1 DT 

- 7 T = ^ - 7 T 

dx T Dt 



^Vg ht f 
rMPcp _ rpgCp 


U h \ 

r\PgCpU L / 
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( 8 ) 


st , = ht f 1 ( h \ 
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The numbers Stg and Sty have the form of Stanton numbers modified by the 
presence of the dimensionless lengths r and Z w , respectively. The use of a 
parameter Sty, containing both fluid and wall properties, is unusual. The ratio 


2s = Z wPw c w 
St w r Pg c p 


which is equal to one-half the ratio of the heat capacity of the wall to the 
heat capacity of the gas, could be used in place of Sty. However, St w has 
been retained since it arises naturally in the development of the equations. 

It is seen that Stg and Sty. completely d_gte^mine the differential equa- 
tions for the dimensionless dependent variables T, Tw.> and u. It is shown in 
reference 2 that the dimensionless constants 


A 




( 10 ) 


T L - T 


( 11 ) 



8 


enter the initial and boundary conditions fpr the dimensionless equations. 
Within the assumptions made thus far Stg, L 0 , and Tp are constant for a 
given problem. The St^ will vary only if c w is allowed to vary. 


Pressurant Mass Ratio 


Defining an ideal pressurant mass 

mj_ = itr 2 (L f - L 0 )p g 

it can be shown that the mass ratio (sometimes called collapse factor) is 
given by 


A A 



(12) 


The mass ratio is, therefore, known when the solution^of equations (4), (5), 
and (6) for the dimensionless temperature variation T(x,^) is known. These 
considerations lead to the following conclusion: With the assumptions stated 

in the analysis, and with the further assumption that the wall specific heat 
is constant, the mass ratio is completely determined by the specification of 
four dimensionless constants St g , Stw, L 0 , and Tp. 


This conslusion is not restricted to any particular liquid, pressurizing 
gas, or tank wall material. The constant L 0 is determined by the initial ul- 
lage ratio, and the constant Tp is determined by the saturation temperature 
and the pressurizing gas temperature. All other characteristics of the problem, 
for example tank wall material, wall thickness, tank radius, density, and specific 
heat of the pressurizing gas and tank pressure, enter only through the constants 
Stg and St^. Within the assumptions of the analysis, therefore, a complete param- 
etric investigation can be done by examining the effects on the mass ratio of 
variations in Stg, St^, L Q , and Tp. 

For hydrogen problems, however, the assumption that c w is constant is not 
very good. If this assumption is dropped, the preceeding conclusion no longer 
holds. The specific heat c w then varies with temperature T w and the form of 
of the variation may change from one wall material to another. This leads to 
the following conclusion: With the assumptions stated in the analysis, and con- 

fining attention to a single wall material, the mass ratio is completely deter-^ 
mined by the specification of four dimensionless constants Stg, St w , L 0 , and Tp 
and the inlet gas temperature T g . 
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Effect of Parameters 

/S 

^ The parametric investigation is then continued as follows. Values of L 0 , 
Tl, and T g are fixed and computer solutions of equations (5), (6), and (7) 
for a wide range of values of St g and Stv are obtained. From these solu- 
tions (in particular, the temperature distributions) the mass ratios are com- 
puted. The results of these calculations are shown in figure 3. 

For fjxed values L 0 = 0.0526 (corresponding to an initial ullage of 5 
percent), Tl = 0.074, and T g = 500° R, figure 3 enables the prediction of 
pressurant mass ratio (collapse factor) for a wide range of design conditions, 
within the assumptions of the analysis. 

A A 

The effect , of the arbitrarily chosen values of Tg, Tl, and L 0 is exam- 
ined next. Representative curves (St g =5.0 and St w = 2.5) are taken from 
figure 3. With these curves fo£ comparison the value of T g is changed to 
300° and 700° R with Tp and L 0 held at their original values. Again mass 
ratios are obtained from computer solutions and the results are compared with 
the original results (fig. 3) for T = 500° R. Figure 4 gives an indication 
of the effect of Tg on the mass ratio. The effect is large only for large 
values of Stg. 

/S /N 

In a similar manner the effect of Tl is found^by holding Tg and L 0 
fixed at the values used for figure 3 and changing Tp to 0.12 (’i'p = 0.12 
corresponds to Tp = 60° R and Tl = 0.074 corresponds to Tp = 37° R) . The 
results are shown in figure 5. The effect on the mass ratio is small. 

A 

The dimensionless initial ullage height L Q is treated similarly, chang- 
ing it from the value 0.0526 (corresponding to an initial ullage volume of 5 
percent) used in figure 3 to the value 0.25 (corresponding to an initial ullage 
volume of 20 percent). As shown in figure 6 the initial ullage effect is small 
for values of initial ullage up to 20 percent. 

Figures 4 to 6 indicate that the reference Stanton number map (fig. 3) 
has a wider range of validity than was first evident. In particular, the uge 
does not appear to be restricted to the particular values of Tg, Tl, and L 0 
that were used to obtain figure 3. This conclusion will be checked against 
experimental data in a later section. 


Effect of Assumptions 

It is possible to examine, in a similar manner, the effect of some of the 
assumptions entering the analysis. Figure 7 shows results obtained using a 
variable gas specific heat. The difference is negligible. Figure 8 shows the 
relatively large effect, on the other hand, of choosing wall specific heat to 
be constant. It was this latter result that led to the inclusion of varying 
wall specific heat in determining the reference Stanton number map. It is in- 
teresting that changing the wall material from stainless steel to aluminum has 
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little effect on the mass ratio (fig. 9). The reason for this is given in 
reference 2 . 

It is shown in reference 2 that the choice of initial values of gas and 
wall temperatures affects the mass ratio little for initial ullages up to 
20 percent. The effects of initial transients in outflow rate and inlet gas 
temperature are shown in that report to he small. Transient pressure effects 
are more important. 


Comparison with Experiment 

An analysis of the tank pressurization problem has indicated that the 
primary parameters affecting the mass required to pressurize a cylindrical 
tank during outflow can be combined into two dimensionless groups having the 
form of modified Stanton numbers , one associated with the gas and one with the 
tank wall. This enables approximate values of mass ratio (collapse factor) to 
be determined from a single figure for a large range of design variables. To 
test this conclusion the experimental data used previously in the paper will 
be used again. 

In the case of the Lewis experiments and the Lockheed-Georgia experiments 
described before, the experimental average values of heat-transfer coefficient 
are available. Using these values of h the gas and wall Stanton numbers can 
be determined for each set of data. Using these Stanton numbers and figure 3, 
an estimated value of mass ratio can be obtained. 

Values of mass ratio determined in this way for the Lewis experiments are 
shown in table IV. One of the Lockheed-Georgia experiments was omitted since 
it contained helium in the initial ullage space and was subsequently pressurized 

TABLE IV. - COMPARISON OF EXPERIMENTAL VALUES OF MASS RATIO WITH VALUES 


DETERMINED FROM THE REFERENCE STANTON NUMBER MAP 


Example 

Experi - | 

Mass ratio 

Percent! Example 1 

Experi- 

j Mass ratio ! 

-Percent 


mental 

determined 

differ-' 

t 

mental 

i determined 

differ- 


■mass ratio i 

from 

ence 

! 

mass ratio! from 

ence 



Stanton number 




Stanton number 




map 


1 


map 


, 

Lewis data 



Lockheed- 

-Georgia data 


1 

2.58 

2.76 

-6.5 

1 

1.72 

! 1.54 ! 

1 11.7 

2 

1.77 1 

1 1.84 

-3.8 

2 

2.14 

2.16 

-0.9 

3 

3.09 

3.31 

-6.7 

3 

1.79 

1.61 

11.2 

4 

3.20 

2.25 

-2.2 

4 

1.71 

; 1.68 

1.8 

5 

1 .47 

1.35 

8.9 

5 

1.81 

1.69 

7.1 

6 

2.37 

' 2.56 

-7.4 

7 

1.80 

1.71 

5.3 

7 

! 2.73 

2.86 

-4.5 

8 

1.75 

1.59 

10.0 

8 

1.86 

1.93 

■ -3.6 

9 

1.83 

1.69 

8.3 

9 

1.38 

1.25 

• 10.4 





10 

3.92 

4.25 

-7.8 






| i 
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with hydrogen, a situation not covered by figure 3. Values of mass ratio 
determined from figure 3 for the other eight Lockheed-Georgia experiments are 
shown in table IV. Also shown are the experimental values of mass ratio for all 
the experiments. In the case of the Lewis data the actual experimental values 
are modified to eliminate the heat sink effect that is not accounted for in 
figure 3. 

The percent difference between calculated and experimental values of mass 
ratio is also shown in table IV. The average value of the absolute differences 
for the Lewis data is about 6 percent. For the Lockheed-Georgia data the average 
is about 7 percent. These results bear out the implications of the parametric 
analy si s . 

It should be remembered, however, that to compute Stanton numbers for 
design purposes a value of heat -transfer coefficient h must be estimated. In 
reference 2 a simple method of estimating h from a free convection formula 
is examined. For the experiments considered here such a simple method appears 
to be adequate. Its general use, however, is open to serious question and the 
determination of heat-transfer coefficient for arbitrary conditions remains an 
unsettled question. 


APPENDIX - SYMBOLS 

C effective perimeter of internal hardware 
Cp specific heat of gas 

c w specific heat of tank wall 

h heat -transfer coefficient 

J mechanical equivalent of heat 
Lp ullage height at time t = tp 

A 

Lp dimensionless ullage height, Lp/(Lf - L 0 ) 

L 0 ullage height at time t = 0 

A 

L q dimensionless initial ullage height, hj - L c ) 
l height of ullage (see fig. l) 

Z w thickness of tank wall 

✓S . 

7 W dimensionless thickness of tank wall, Z w /(Lf - L c ) 

M molecular weight 
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m mass of pressurant gas added during outflow 

m i mass of pressurant gas required assuming no heat transfer 

P pressure in tank 


'll 


*0 

R 


r 

A 

r 


St g 


Stw 

T 


T 



t 


% 

t f 

At 

u 


heat flow rate to gas from internal hardware 
heat flow rate to tank wall from outside 
universal gas constant 
radius of tank 

dimensionless radius of tank, r/(L^. - L ) 

modified gas Stanton number, i — 

r,.Pg c pU L 

modified wall Stanton number, ^ — 

^w Pw c w u L 

gas temperature 

dimensionless gas temperature, T/Tg 
gas temperature at tank inlet 

gas temperature at liquid interface 

dimensionless temperature, Tp/Tg 

temperature of tank wall 

dimensionless temperature, T w /Tg 

time 

dimensionless time, t/t^ 
time at end of outflow 

time increment for finite difference equations 
velocity of gas 


u dimensionless gas velocity, u/u^ 
u-^ velocity of gas at liquid interface 
x space coordinate in direction of tank axis 

x dimensionless space coordinate, x/(Lf - p ) 


13 


Ax space increment for finite difference solution 

Z compressibility factor 

Pg density of gas 

p w density of tank wall 
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Figure 2. - Comparison of calculated and experimental gas and wall temperatures at end of outflow. 
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Figure 1. - Schematic drawing of cylindrical tank. 
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Figure 2. - Concluded. Comparison of calculated and experimental gas and wall temperatures at end of outflow. 
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Figure 3. - Stanton number map showing values of mass ratio for range of gas and wall Stanton numbers. Initial 
ullage ratio, 0.05; dimensionless interface temperature, 0.074. 
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Figure 4 - Effect on mass ratio of changing the inlet gas temperature. Dimensionless interface temperature, 0.074 
dimensionless initial ullage height, 0.0526. 




(b) Modified wall Stanton number, 2.5. 

Figure 5. - Effect on mass ratio of changing the dimensionless interface temperature. Inlet gas temperature, 500° R; 
dimensionless initial ullage height, 0.0526. 



.01 .02 . 04 . 06.08 .1 .2 . 4 . 6 . 8 1 2 4 6 8 10 20 40 60 80 100 

Wall Stanton number, 


(a) Modified gas Stanton number, 5. 



<bl Modified wall Stanton number, 2 . 5 . 


Figure 6 . - Effect on mass ratio of changing the dimensionless initial ullage height. Inlet gas temperature, 500° R; 
dimensionless interface temperature, 0 . 074 . 
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(b) Modified wall Stanton number, 2.5. 

Figure 7. - Effect on mass ratio of allowing gas specific heat to vary with temperature. Inlet gas temperature, 500° R; 
dimensionless interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 
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(b) Modified wall Stanton number, 2.5. 


Figure 8. - Effect on mass ratio of holding wall specific heat constant Inlet gas temperature, 500° R; dimension- 
less interface temperature, 0.074,- dimensionless initial ullage height, 0.0526. 
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(b) Modified wall Stanton number, 2.5. 

Figure 9. - Effect on mass ratio of changing tank wall temperature from stainless steel to aluminum. Dimensionless 
interface temperature, 0.074; dimensionless initial ullage height, 0.0526. 
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